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Abstract. The stability of magnetized strange quark matter (MSQM) is investigated 
within the phenomenological MIT bag model, taking into account the variation of the 
relevant input parameters, namely, the strange quark mass, baryon density, magnetic 
field and bag parameter. We obtain that the energy per baryon decreases as the 
magnetic field increases, and its minimum value at vanishing pressure is lower than the 
value found for SQM. This implies that MSQM is more stable than non-magnetized 
SQM. Furthermore, the stability window of MSQM is found to be wider than the 
corresponding one of SQM. The mass-radius relation for magnetized strange quark 
stars is also derived in this framework. 
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1. Introduction 

Strange quark matter (SQM), in contrast to ordinary nuclei (where quarks are confined 
into colorless nucleons), is thought to be a form of macroscopic matter which contains a 
large quantity of deconfined quarks (m, d and s) in /5-equilibrium, with electric and color 
charge neutrality. Theoretically, Bodmer [1] and Witten [2] proved that SQM could be 
the more stable phase of nuclear matter. The reason behind this crucial conclusion is 
the presence of the quark s with a mass lower than the Fermi energy, which leads to a 
decreasing of the total energy per baryon E/A, if some of the u and d quarks on the 
Fermi surface are converted into s quarks via weak interactions. This is, in other words, 
a consequence of the Pauli exclusion principle. 

If the hypothesis of the existence of SQM is true, this form of matter could have 
interesting astrophysical consequences. For instance, SQM could be realized in the 
inner core of compact objects, such as neutron stars, or even form the hypothetical 
strange stars. It may be possible as well that small quarks droplets could be formed in 
relativistic heavy ion collider (RHIC) experiments [3]. As SQM is more stable and bound 
at finite density, one expects configurations of quark stars with macroscopic properties 
quite different from neutron stars [H |5]. These objects would be self-bound and their 
masses would scale with the radius as M ~ i?^, in contrast to neutron stars which 
have masses that decrease with increasing radius (M ~ R~^) and are bound by gravity. 
Thus, they could explain a set of peculiar astrophysical observations which are awaiting 
for theoretical explanations P, [3, El [9l [IDj ■ 

SQM has been widely studied using the MIT bag model, where the non-interacting 
quark gas at zero temperature is considered to be confined into a phenomenological 
bag, -Bbag, which mimics the strong interaction. Such studies have revealed that even 
heavy strangelets might be stable, with the corresponding E/A lower than that of the 
iron nuclei, E/A{^^Fe) ~ 930 MeV [H]. Being a phenomenological model, the MIT 
bag model has some limitations that are well known. The confinement is put by hand, 
and the model is unable to reproduce the chiral symmetry breaking at zero density. 
Nevertheless, the model is quite satisfactory in describing quark matter at finite density. 

Focusing on a real astrophysical scenario, one should also note the important role 
played by the strong magnetic fields expected in this context. Pulsars, magnetars, 
neutron stars, the emission of intense sources of X-rays could be associated to sources 
with intense magnetic fields around 10^^ — 10^^ G or even higher fields [I2l [13]. We 
remark that the magnetic field intensity may vary significantly from the surface to the 
center of the source. Although there is no direct observational evidence for the magnetic 
field strength in the inner core of stars, theoretical estimates indicate that fields as high 
as 10^^ G could be allowed. Therefore, SQM in the presence of strong magnetic fields 
is a subject that it is worth addressing. In a recent work [Tl], we investigated the SQM 
properties in the presence of a strong magnetic field, taking into account (3- equilibrium 
and the anomalous magnetic moments (AMM) of the quarks. It was found that the 
stability of the system requires a magnetic field value B < 10^^ G. This is in contrast 
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to the bound B < 10^^ G, obtained when AMM are not included [THj. 

The scope of the present paper is to estabhsh a bridge between the microscopic 
properties of magnetized SQM (MSQM), given by the equation of state (EoS) of the 
system, and the macroscopic observables of strange quark stars (SQS), derived from 
this kind of magnetized matter. In the framework of the MIT bag model, we study the 
behavior of the system with the variation of its parameters, namely, the bag parameter 
-Bbag; the strange quark mass m^, the baryon density and the magnetic field B. The 
possible mass-radius configurations of magnetized strange quark stars (MSQS) are then 
obtained by solving the Tolman-Oppenheimer-Volkoff (TOV) equation. As it turns out, 
stable stars with smaller radii are allowed due to the compactness of matter and the 
presence of a strong magnetic field, since the energy per baryon at vanishing pressure is 
lower in this case. 

The paper is organized as follows. In section [2] we briefly review the thermodynamic 
properties of SQM in the presence of a strong magnetic field. In section [3] the stability 
windows of SQM in the presence of a magnetic field are obtained varying the relevant 
input parameters of the model. Section H] is devoted to the study of the mass-radius 
relation for MSQS by numerically solving the TOV equations. A comparison with some 
of the available observational data [16] is also presented. Finally, our conclusions are 
given in section [5l 



2. SQM in the presence of a strong magnetic field 

In a previous work [S], the thermodynamical properties of MSQM were studied within 
the MIT bag model and including the AMM of quarks. Assuming that quarks freely 
move inside the bag, the total thermodynamical potential is given by 

^total = ^vac + ^ , (1) 

where Qyac = B^agV is the QCD vacuum energy density, V is the volume and 

n = -T\nZ = ^{niyV + ni^sS + ni^cC), i = e,u,d,s. (2) 

i 

The quantities flv,s,c are the volume, surface and curvature contributions, respectively. 
Terms associated to the surface and curvature play an important role in processes related 
to phase transitions, like bubble nucleation in dense matter. Since for the purpose of 
studying the thermodynamical properties and EoS of SQM in /3-equilibrium only the 
bulk quark matter is relevant, these two terms can be safely neglected. The explicit 
form of the volume contribution is 

a,y(T,/..) = -^I^E / dhHl + e^^>^^-% (3) 

where n denotes the Landau level, /Xj is the chemical potential of i-particle; d^ = 1 
and du^^s = 3 are degeneracy factors. The quantity corresponds to the energy of the 
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constituents, 



,2 _ T 2 



B 



kl + m^J—{2n + l-r^) + l~r^y^B] , (4) 



where -Bf = mf/\ei\ is the critical magnetic field, = \Qi\/mi accounts for the AMM 
and ?7 = ±1 correspond to the orientations of the particle magnetic moment, parallel or 
antiparallel to the magnetic field. 

For a degenerate MSQM, the energy density, pressures and number density are 
given by the expressions [H] 

^ = 5 E E (-^P%^ + hf ' In "-^) , (5) 

i n \ * / 

^1 = 5 E E [-^p% - " 1- ' (6) 
= ^ E -^^ E i^^hf 1- ^41^ 1 ' (7) 



i * n 

where the different expressions for the parallel pressure P\\ and the transverse pressure 
reflect the anisotropy of pressures due to the magnetic field. In equations ©-(jH]) we 
have deflned the dimensionless quantities 




Xi = fii/nii, hi = \l — {2n + 1 - 1]) + 1 - i]yiB , 

Ph-^^', ,;^^^^JBi2=^^-,„.B, (9) 

where Xj is the dimensionless chemical potential, pp^i corresponds to the modified 
Fermi momentum due to the magnetic field and hi corresponds to the magnetic 
mass. The sum over the Landau levels n is up to nl^^^ given by the expression 
'^max = ^ [((^i + VViBy — 1) B^ / {2B)], whcrc I[z] denotes the integer part of z. 

Let us remark that the quantities defined in (I5l)-(l9]) contain the contribution of 
Landau diamagnetism (given by the quantization of Landau levels) as well as the Pauli 
paramagnetism (due to the presence of AMM for quarks) [Hj. Since the inclusion of 
AMM will not change our main conclusions (it just places a more restrictive upper 
bound on the magnetic field), in what follows we shall neglect their contribution and 
consider only the effect of Landau diamagnetism. This can be easily done by taking 
j/j = in all the expressions. 
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Figure 1. Energy per baryon as a function of the baryon density for B = and 
B = 5 X 10^^ G, assuming the stability condition of a vanishing pressure. 



3. Stability window for magnetized SQM 

In this section we study the stabihty window of SQM in a strong magnetic field as a 
function of the input parameters of the model: the baryon density hb, the magnetic 
field B, the bag parameter Bf,ag and the quark mass mg. In the context of the MIT bag 
model, and in the absence of a magnetic field, the stability condition for SQM requires 
the study of the equations 



under the condition P = 0, where P and e are the total pressure and energy, respectively. 
For SQM with massless quarks one obtains the well-known EoS P = {e—4Bbag)/3, which 
at vanishing pressure then leads to the stability condition e = 4:Bhag. 

In the presence of a strong magnetic field, the anisotropy in the pressures implies 
P± < P||. Thus, the stability condition for strong fields changes from P = to P^ = 
or, equivalent ly. 



As it turns out from our numerical analysis (see figures below), in this case the total 
energy e given in ffTOl) is always lower than ABfyag and, therefore, MSQM is more stable 
than non-magnetized SQM. 

In order to obtain the EoS of MSQM in the stellar scenario under consideration (i.e. 
the inner core of neutron stars or quark stars) we should also take into account a set of 
equilibrium conditions: /3-equilibrium, baryon number conservation and electric charge 
neutrality. Lepton number is not conserved because neutrinos are assumed to enter and 




(10) 
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Figure 2. Energy per baryon versus pressure for S = and B — 5 x 10^^ G. 
We take Sbag — 75 MeV/inx^. The horizontal dot-dashed hne corresponds to 
E/A{^^Fe) ~ 930 MeV. 




Figure 3. EoS: Pressure versus energy for B — and B = 5 x 10^® G. For comparison, 
the SQM EoS with massless quarks, i.e. P = (e — 4_Bbag)/3, is also shown. 

leave the system freely (/Xj^ = 0). The equilibrium conditions are then determined by 
the following relations: 

2Nu -Nd-N,-?,Ne = Q, (12) 

The system of equations (fTT|) - (fT2|) involves all the parameters of the MSQM model. The 
system is completely determined: there are five equations to find the quark and electron 
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Table 1. Bag parameter and baryon density in the presence of a strong magnetic field. 
We assume = 150 MeV and E/A = 930 MeV. 





5bag (MeV/fm3) 


HB/no 





75 


2.18 


5 X lO^s 


80 


2.37 



chemical potentials and the bag parameter Bbag- 

In figured] we present a comparison of the energy per baryon E/A (or, equivalently, 
e/ub) versus the number density n^/rio ('^o — 0.16 fm~^) in the absence of a magnetic 
field and for a magnetic field value of 5 x 10^^ G, both at vanishing pressure. We assume 
rriu = TUd = 5 MeV and nig = 150 MeV. As can be seen from the figure, E/A is always 
lower in the presence of a magnetic field. Figure [2] shows the behavior of the energy per 
baryon E/A with the pressure P, for B = and B = 5x 10^^ G. The bag parameter has 
been chosen to be 5bag = 75 MeV/fm^. One can notice that the point of zero pressure 
for MSQM is reached for an energy density value lower than that of SQM. Consequently, 
matter is indeed more stable and more bound when a magnetic field is present. At the 
zero-pressure point the corresponding baryon density is ub — 2. 18 no for i? = and 
ub ^ 2.08^0 for 5 = 5 X 10^^ G. 

Aiming at obtaining the MSQS observables (e.g. the star mass and radius), next we 
determine the EoS of MSQM. Since in the present case we do not have a simple algebraic 
relation between the pressure P and the energy density e, such a relation must be found 
numerically. In figure [3] we show the EoS for SQM (i.e. when B = 0) and for MSQM for 
a magnetic field around 5 x 10^^ G. For comparison, we also include the EoS of SQM 
neglecting all the quark masses: P = {e — 4i?bag)/3. Although the changes in the EoS 
are not too significant, they can alter the macroscopic observables of stars formed by 
MSQM, as we shall see in the next section. 

Let us now study the stability window for strange quark matter and, in particular, 
the regions in the (m^, ?T,B)-plane where SQM is stable. In order to investigate how the 
magnetic field affects this window, we fix the magnetic field to -B = 5 x 10^^ G and study 
the contours of 5bag and E/A. The results are then compared with the ones obtained 
for SQM at -B = 0. In figure H] we present the contours of constant bag parameter and 
energy per baryon for SQM. The corresponding contours for MSQM are depicted in 
figure [5l 

As can be seen from the figures, the magnetic field tends to shift the stability 
window of SQM towards higher values of the baryon density (cf. table [I]). Indeed, 
for a given strange quark mass, the constant lines of i?bag and E/A in figure H] are 
displaced to a higher value of ub/uo, when compared with the corresponding one of 
figure [5l While for SQM the allowed range for the baryon density is 1.8 ^ '^b/'^o ^ 2.4 
for 50 < < 300 MeV and Ebag > 57 MeV/fm^, MSQM allows densities in the 
range 1.85 < riB/riQ < 2.6 for 50 < nig < 240 MeV. The approximate lower bound 
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of 57 MeV/fm^ for the bag parameter is the hmit imposed by requiring instabihty of 
two- flavor quark matter (deconfined gas of u and d quarks) [HI [17]. We also notice that 
the allowed interval for the bag parameter is affected by the magnetic field. Below the 
energy contour of 930 MeV, our EoS for MSQM corresponds to an E/A at P = lower 
than that of ^^Fe and it can therefore be considered absolutely stable. 
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Figure 4. The stability windows of SQM in the (ms, 7iB)-plane. The contours of 
constant .Bbag (dashed hnes) and E/A (soUd hnes) are shown. 




Figure 5. The stabihty windows of MSQM in the (m^, ris) -plane for S — 5 x 10^^ G. 
The contours of constant i?bag (dashed lines) and E/A (solid lines) are shown. 
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4. SQS mass-radius relation 

In this section we study the equilibrium configuration of magnetized strange quark stars 
described by the EoS of MSQM obtained in the previous section. As is well known, the 
most important macroscopic parameters of a star are its radius R and its gravitational 
mass M. The surface redshift Zg is also an important parameter which is related to the 
star mass and radius through the equation 

^ - 1, (13) 



1 2GM 



where G is the gravitational constant. The interest in this parameter is that it is 
observable. 

Configurations of spherical symmetric non-rotating compact stars are obtained by 
the numerical integration of the TOV equations 
dM 



dr 



AirGr'^eir 



dP _ ^ {ejr) + P(r)) (M(r) + AnPjry) ^^^^ 
dr — 2rM{r) ' 

supplemented with the EoS, where P(r) is the pressure and 6{r) is the energy density. 
The radius R and the corresponding mass M of the star are determined by the value 
of r for which the pressure vanishes, P{R) = 0. The EoS fixes the central pressure, 
P(0) = Pc, which together with the condition M(0) = 0, completely determine the 
system of equations (fT^ . In this way, varying continuously the central pressure, one 
obtains a mass-radius relation M{R), which relates masses and radii for a given EoS. 
The stable branches of these curves must satisfy the condition dM/dPc > 0. Other 
solutions are unstable and collapse. 

Let us now study the static stellar configurations for MSQS. We remark that, 
despite the anisotropy of pressures in the presence of the magnetic field, this effect is 
not significant for B < 10^^ G. Therefore, in our astrophysical context, even if the EoS 
is determined by the condition ffTTl) . TOV equations can be solved using a spherical 
metric, which then leads to equations dHl). 

Starting from the numerical EoS of MSQM, we have solved the TOV equations 
to obtain the stable mass-radius configurations of MSQS. As mentioned before, the 
magnetic field enhances the stability of SQM and, thus, it allows the appearance of 
stable stars with masses and radii smaller than for non-magnetized SQM. 

In figureOwe plot the stable configurations of SQM without a magnetic field, as well 
as of MSQM, for two different values of the magnetic field, B = 5 x 10^^ G and 10^^ G. 
The curves are presented for two values of the bag parameter, -Bbag = 57 MeV/fm^ and 
75 MeV/fm^. Our main results are summarized in tables [2] and [3l where the maximum 
mass Mmax and maximum radius -Rmax of MSQS are given. 

For the sake of comparison with observational data, we have also depicted in figure [6] 
the la and 2a contours for the mass and radius of the low mass X-ray binary EXO 1745- 
248, based on the spectroscopic data during thermonuclear bursts combined with a 
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R[km] 

Figure 6. Possible M—R configurations obtained for MSQS in the present model. The 
curves are shown for different values of the magnetic field B and bag parameter -Bbag 
(in units of MeV/fm^). The shaded (blue) areas correspond to la and 2a confidence 
contours for the mass and radius of the X-ray binary EXO 1745-248 [16j . 

Table 2. Maximum mass and the corresponding radius of a MSQS for different values 
of the magnetic field and bag parameter. 



B{G) 


-Bbag — 


57 


-Bbag — 


75 




R (km) 




R (km) 





1.87 


10.35 


1.65 


9.07 


5 X IQi® 


1.84 


10.16 


1.62 


8.91 


1019 


1.76 


9.73 


1.56 


8.62 



Table 3. Maximum radius and the corresponding mass of a MSQS for different values 
of the magnetic field and bag parameter. 



B{G) 


Sbag — 


57 


-Bbag 


= 75 


-Rmax (km) 


M/Mq 


i?max (km) 


M/Mq 





10.80 


1.73 


9.47 


1.49 


5 X 


10.60 


1.69 


9.32 


1.50 


1019 


10.17 


1.58 


9.00 


1.41 



distance measurement to the globular cluster [16]. We observe that the curves predicted 
by the MSQM EoS are in good agreement with the confidence contours. 
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5. Conclusions 

In the present work, we have investigated the stability of magnetized strange quark 
matter within the phenomenological MIT bag model. We have studied the stability 
windows of MSQM taking into account the variation of the strange quark mass, the 
baryon density, the magnetic field and the bag parameter. We found that MSQM 
is indeed more stable than non-magnetized SQM. While for SQM the stability range 
for the baryon density is 1.8 < ub/uo < 2.4 for 50 < < 300 MeV, MSQM 
allows densities in the range 1.85 < ns/riQ < 2.6 for 50 < < 240 MeV and a 
magnetic field value of 5 x 10^^ G. Moreover, the allowed range for the bag parameter 
is 57 < Sbag < 90 MeV/fm3. 

At vanishing pressure and finite density, the derived EoS for MSQM exhibits a 
minimum of the energy per baryon lower than that of ^^Fe. This gives the possibility 
of existence of stable configurations for magnetized strange quark stars. In this 
simple framework for MSQM, we have then derived the mass-radius relation for such 
compact stars and compared the theoretically predicted relations with some of the recent 
observational data. We have concluded that in the presence of a strong magnetic field, 
stable strange stars with smaller masses and radii could exist. 
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